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Using a functional renormalization group method, we obtain the phase diagram of the two-leg lad- 
der system within the Holstein-Hubbard model, which includes both electron-electron and electron- 
phonon interactions. Our renormalization group technique allows us to analyze the problem for both 
weak and strong electron-phonon coupling. We show that, in contrast results from conventional weak 
coupling studies, electron-phonon interactions can dominate electron-electron interactions because 
of retardation effects. 
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The renormalization group (RG) method has provided 
fundamental understanding of the stable phases of mat- 
ter in terms of the basin of attraction of fixed points 
(FP). Quantum phase transitions (QPT) between differ- 
ent fixed points can be described in terms of the stability 
of these FP relative to the RG flow. In the theory of 
metals, Landau's Fermi liquid theory can be described 
either as a stable FP of the renormalization group flow 
under repulsive interactions, or as an unstable FP un- 
der attractive interactions to a superconducting state 
(SC) described by BCS theory Apart from deter- 
mining the stability of a FP, the RG method also pro- 
vides direct information about the energy scales (gaps) 
and length scales (correlation lengths) associated with a 
QPT. Hence, the RG not only provides qualitative in- 
formation, but also quantitative information, about the 
different phases of matter. 
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FIG. 1: Phase diagram of Holstein-Hubbard ladder as a func- 
tion of electron-phonon coupling g^p and doping {t = 1,U = 2, 
Wo = 1)- Zero doping corresponds to half-filling and doping 
of 0.5 corresponds to quarter-filling. 

In practice, however, there are very few examples of 
problems in which the RG can be implemented beyond 
weak coupling 2]. Furthermore, another well-known 
challenge to the RG has been its implementation in prob- 
lems with retardation, such as the electron-phonon (e- 



ph) problem, because of the introduction of multiple, 
frequency-dependent interaction channels that renormal- 
ize in a collective, energy dependent, way. At weak e-ph 
coupling, when retardation plays a small role, the RG 
can be implemented using approximations such as the 
so-called " two-step" RG 0, 01 • Usually, the two-step RG 
leads to an enhancement of the instabilities that exist 
when only electron-electron (e-e) interactions are present, 
but does not describe phases where retardation effects 
due to the e-ph coupling become dominant over e-e ef- 
fects. In fact, we show that new phases, such as charge 
density wave (CDW) and s-wave SC, emerge at strong 
e-ph couphng (see Fig.P). 

Recently, the RG was extended to study the non- 
perturbative physics of the e-ph problem in the large 
N limit (where N = Ep/A, Ep is the Fermi energy 
and A is the cut-off) in the presence of electron-electron 
(e-e) interactions for a two-dimensional circular (three- 
dimensional spherical) Fermi surface 0. In particular, 
that work established that the stable FP of this prob- 
lem is described by the Eliashberg-Migdal theory, which 
contains BCS as its weak-coupling limit. Since it has few 
scattering channels P| , the circular (spherical) Fermi sur- 
face, which is a consequence of the plane- wave character 
of the electronic wavef unctions, is a particularly simple 
problem. This is not the case for Fermi surfaces that 
derive from localized electronic orbitals that can be ei- 
ther unstable to exotic SC (such as p- and d-wave) or 
to density wave (charge and spin) phases. Arguably, the 
most famous example of this phenomenon are the one- 
dimensional (ID) conductors where the nesting associ- 
ated with Fermi points leads to a richness of phases that 
can be described^ by g-ology, bosonization 0,13], and the 
functional RG 8]. However, the ID problem is very par- 
ticular because of the restricted phase space for interac- 
tion, and even in ID, the treatment of the e-ph interac- 
tion is problematic jsj ■ 

The recently developed e-ph coupled functional renor- 
malization group (FRG) y includes the merits of the 
FRG for the pure e-e interacting system P, 0, 0| and 
goes beyond the two-step renormalization group for the 
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e-ph interacting system 0, 0, 0, by allowing the sys- 
tematic study of the full frequency dependence of cou- 
plings. In this paper we use this newly extended FRG to 
study the two-leg ladder Holstein-Hubbard model that 
describes a quasi-lD system in the presence of e-e and e- 
ph interactions. The ladder system, consisting of two ID 
chains coupled by inter-chain hopping, can be considered 
as an intermediate between ID and 2D systems, since 
it shares properties with both: it has a limited number 
of scattering channels, as in ID, but also can show SC 
with 2D characteristics, such as d-wave order parameter 
[l^ H^ . As a result, ladders have been considered by 
many as a theoretical testbed for the understanding of 
high-Tc superconductivity (111] . 

Furthermore, the interest in the e-ph problem in high- 
Tc has been revitalized by experimental evidence that 
the e-ph coupling may play a critical role in these systems 
0,0]. Apart from the cuprates, ladder systems have 
also attracted much attention recently because, from the 
technical point of view, they are amenable to field theo- 
retical calculations and to reliable numerical simulations, 
in contrast to truly 2D systems. Studying the inter- 
play between the e-e and e-ph interactions in a ladder 
is also important for the understanding the physics of 
some low-dimensional materials, such as molecular crys- 
tals and charge transfer solids |18j . 

The Holstein-Hubbard model (HHM) is described by 
the following Hamiltonian (we use units such that h ~ 
l = kB): 

H = ~t {c\ „Cj^„ + H.c.) + UY^ ni^^m^i 

+9cpY{a\ + ai)ni^„ + ujQ^alai , (1) 

where c\ ^ {ci,a) creates (annihilates) an electron at site 

i with spin a {nia is the electron number operator), a\ 
{tti) creates (annihilates) an optical phonon at site i with 
energy wq, t\s the hopping energy of the electrons along 
legs and rungs of the ladder, U is the on-site e-e interac- 
tion, and gop is the e-ph coupling energy. In the absence 




FIG. 2: Energy dispersion for free electrons in a 2-leg ladder. 



of e-e (L'^ = 0) and e-ph (^ep = 0) interactions, the energy 
spectrum of the electrons has the form shown in Fig. |2] 
In what follows we use units such that t = 1. 

One can reformulate the problem in terms of path in- 
tegrals, and since the phonon fields enter quadratically 
in the action, they can be integrated out exactly , gen- 
erating an effective action 5" for the electrons that can be 
written as: 

S ^ [ ^Umu; - ek)Mk) (2) 

J{k}.,cr,a' 

where ip is the electron field, and the coupling functions 
for electron-electron interactions in the momentum space 
are given by: 

g{ki,k2,h,k4) = U ~ 2glpUJa/[ujl + {uji - ^4)^], (3) 

where k — {k,uj). Each coupling is a function of the 
in-coming and out-going momenta and frequencies, sub- 
jected to the conservation of momenta, ki+k^ = fca + k^ 
and energy oji + uj2 = oJz + loa- In a purely e-e interact- 
ing system, the coupling functions are always indepen- 
dent of the frequency. For a generic system with e-ph 
interactions, the couplings become frequency dependent, 
i.e., retarded. We choose 5(A^l_, fe, fcs, ^4) to represent the 
scattering between electrons with opposite spins. The 
coupling functions for scattering between electrons with 
parallel spins are determined by those with opposite spin 
due to the SU(2) spin symmetry. 

Following ref. Q , our only assumption is that the full 
bandwidth of the two-leg ladder system is larger than any 
of the interaction terms: i.e., that 6t ^ J7, gep, wq, leaving 
the relative values of U, and g^^, free. The strong cou- 
pling limit of the electron-phonon problem occurs when 
A = 2N{0)g'^p/ujQ 3> 1. This regime is accessible in our 
functional RG approach. Utilizing the time-reversal sym- 
metry, exchange symmetry, and reflection symmetry, the 
number of independent couplings in the Fermi surface for 
the 2-leg ladder at incommensurate fillings is 12 (each 
coupling is a function of three momenta and three fre- 
quencies, the frequency dependences are suppressed for 
clarity). There are 6 intra-band couplings: 

9i — J kf;, fc^, fc^), g^ — g{kf^, k^, k^, k^), 

91 — J k^, fc^, fc^), gi ~ g(k^i k^, k^, k^), 

92 — 9{kLTk^,k^,kf^), g2 — g{k^,k^,k^,k^); (4) 

and 6 inter-band couplings: 

9i — ^(^L I I J kf), g^ — gikf;, k^, k^, fc^), 
9i = sl^L J k^,ki^ , fc^), g2 = g{k^, k^,k^, fcf ), 
9i = 9{kL^kfi,k]^,kfi), 92 = 9{kL jk^jk'^jkf). (5) 

For pure e-e interactions the Fermi velocities are renor- 
malized by g^ and gf , and the renormalizations of scat- 
tering between particles on the same side of the Brillouin 
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zone at different bands, and g^, are absent in the 
leading order. This is not necessarily valid for couplings 
with frequency dependence. Moreover, these couplings 
generate self-energy corrections, which have been shown 
to be iniportant for studying the system with retardation 
effects |5j , leading to the renormalization of Fermi veloc- 
ities and quasiparticle lifetime. Thus we keep these cou- 
plings even though they are not renormalized in the non- 
retarded system . In leading order, our renormalization 
group equations for the coupling function, gAikiT ks): 
and for the self-energy, EA(fc), are given by: 

- j dppt^GA {p)Ga ik)]9A ikiM,k)9Aip,k,k3) 

-J dpptlGA {p)Ga (£i)]ffA (£,fc2,gi )gA(fci ,gi ,fc3) 

-jdpptlGA {p)Ga (£2)] {ki,p,q2)gA {qiMM) 
+ 9A{p^,h,q2)gA{q2,k2,k3)+ ,(6) 

dA^A{k) = -J dpdA[GA{p)][2gA{p,k,k)-gA{k,pA)], (7) 

where k = k2_+h^-p^,qi=P + k3-ki,q2^P + k3 
J dp = J dp^^ l/(27r/3), and Ga is the self-energy cor- 
rected propagator at cutoff A. The initial condition for 
the coupling functions is given by l(2Jl. The self-energies 
are equal to zero at A = Aq. 

For a non-retarded system the low energy instability 
can be obtained from the RG flow of the couplings, and 
different phases can be identified by the FP correspond- 
ing to the different spin and charge modes. This proce- 
dure may not be appropriate for system with electron- 
phonon coupling. Numerical evidence indicates that the 
half-filled Holstcin-Hubbard model violates the Luttinger 
relations 19] . Therefore, instead of following the flows 
of the couplings, we explicitly construct the flows of the 
susceptibilities of different order parameters to obtain a 
unbiased answer within the FRG. The SC susceptibility, 
X^{k,uj), is defined as: 

xi(0,0)=y'i?(l,2)/^(pi)/^(p2)(cp,aC-p,,Tclp,,T4.,i)'(8) 

where 6 — s,d for s- and d-wave SC, respectively [20| . 
The RG equations are: 

Sax! (0,0) = / dpdA[GA{p)GA{-p)]{Zi{p)f, (9) 

9AZi(p)=Yrfp'9A[GA(pOGA(-pO]^V).9A(y,-p:"P,p)ao^ 

The function Zj^[p) is the effective vertex in the defini- 
tion for the susceptibility Xa- Its initial condition, at 
A = Aq, is 1 for s-wave pairing and cos(kx) — cos{ky) 



for d-wave pairing. The RG equations for susceptibilities 
are solved with initial condition Xa=Ao('^'0) ~ "^^^ 
dominant instability in the ground state is given by the 
most divergent susceptibility by solving the RG equa- 
tions numerically. In the calculations reported here, we 
have discretized the frequency axis into 9 points and have 
checked that our results are essentially unchanged when 
additional discretization points are added. 

Without e-ph coupling the d-wave pairing is the generic 
dominant instability for to 0.5 doping (half- to quarter- 
filled), except in a narrow region in the parameter space 
where all the charge and spin degrees of freedom are gap- 
less. The large region of d-wave pairing in the ladder 
comes from the anisotropy between the transverse and 
the parallel directions, which leads the pair scattering to 
develop a directional dependence when the high energy 
degrees of freedom are eliminated. For the Holstein cou- 
pling, the phonon-mediated attractive retarded interac- 
tion is independent of the momentum transfer and is thus 
isotropic. The attractive interaction is relevant in gener- 
ating pair-scattering, since it leads the particle-particle 
ladder diagrams to acquire logarithmic divergences. The 
d-wave pairing from repulsive interaction does not come 
from this divergence, and thus the attractive interac- 
tion is very effective in leading the pairing instability, 
and more importantly, this pairing has s-wave symmetry. 
Therefore, the e-ph coupling tends to enhance the s-wave 
component in the pairing scattering matrix, and once this 
coupling exceeds a critical value, the s-wave susceptibil- 
ity overcomes the d-wave, causing a transition from d- 
to s-wave SC pairing symmetry. This result can only be 
seen in strong coupling and, hence, is not accessible to 
the two-step RG [j| . A study of the fully retarded prob- 
lem for a 2D circular Fermi surface has similarly shown 
how the two-step RG breaks down in the strong coupling 
(A > 1) limit 0. 
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FIG. 3: The flows of s-wave SC (solid lines), and d-wave SC 
(dotted lines) susceptibilities as a function of I = ln(Ao/A) 
for U = 2, uo = 1, and doping of 0.2. 

In Fig. 121 we show a representative RG flow in which 
we observe the change in the dominant susceptibility 
for increasing e-ph coupling g^p. Although the isotropic 
phonons tend to break the d-wave pairs, as long as 
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the e-ph coupling is small compared to the e-e interac- 
tions, the d-wave pairing remains the dominant instabil- 
ity. However, the s-wave overcomes the d-wave pairing 
when the effective interaction in the anti-adiabatic limit, 
C^effectivo ^ U - 2gl^/uJo, IS near zero. As the doping 
increases, the critical values for the pairing transition 
remain close to the points where the anti-adiabatic ef- 
fective interaction vanishes. Notice that for a ladder, the 
d-wave SC region increases with doping and is rather sta- 
ble against s-wave near quarter filling, because the Fermi 
surface becomes more anisotropic with doping. This re- 
sult should be contrasted with the 2D case where d-wave 
superconductivity seems to be suppressed at larger dop- 
ing because the Fermi surface becomes more isotropic. 
Hence, our results show that one has to be cautious when 
comparing results for ladders with the true 2D problem. 
For doping close to the quarter filling, we find noticeable 
changes in the critical values, and a substantially stronger 
electron-phonon interaction is needed to drive the s-wave 
pairing over the d-wave pairing. This can be understood 
from the fact that the difference between the Fermi ve- 
locities of the anti-bonding and bonding bands becomes 
larger as the system is doped away from half-filling (see 
Fig. 121) . From the perspective of the Fermi surface renor- 
malization group, this difference leads to different con- 
tributions to pair scattering at fcy = and at ky — tt, 
which in turn lead to a phase difference in the pairing 
matrix elements which is manifested as a d-wave pair- 
ing channel. This is the reason that the d-wave channel 
becomes stronger as the doping gets closer to quarter fill- 
ing. Nevertheless, the s-wave remains still the dominant 
instability for sufficiently strong e-ph coupling. 

Exactly at quarter filling, the Fermi surface in the anti- 
bonding band shrinks to a point, which leads to a van 
Hove singularity at the center of the anti-bonding band. 
Also, at this filling the distance between the Fermi points 
in the bonding band is commensurate, which results in 
the presence of an extra Umklapp scattering in the bond- 
ing band, — (kf, kf,kj^,k^). We do not find that the 
Umklapp scattering is sufficient to drive the system to 
spin density wave (SDW ) ph ase. In contrast to the case 
without phonons A = [lj| we do not find a Luttinger 
liquid phase (C2S2). The d-wave pairing is still the dom- 
inant instability for small e-ph coupling, and the s-wave 
again overcomes the d-wave at strong coupling. Below 
quarter filling, the anti-bonding band is empty. Only two 
Fermi points remain on the bonding band. The system 
becomes essentially ID as the effect of the anti-bonding 
band is then irrelevant. 

We have also investigated the problem exactly at half- 
filling. In this particular commensurate doping there is 
a direct competition between a SDW and a CDW phase. 
In the absence of e-ph interaction the SDW phase (with 
a charge gap but no spin gap) is very robust and is not 
destabilized at weak coupling. At intermediate coupling 
(A « 1), however, the SDW phase becomes unstable to 



a CDW phase with a charge gap. This result cannot be 
obtained in the weak coupling methods discussed previ- 
ously. 

The final results of our study are summarized in the 
phase diagram of Fig. We find that the e-ph coupling 
at intermediate and strong coupling leads to substan- 
tial qualitative differences from previous results, includ- 
ing the appearance of new phases: CDW and s-wave SC. 
Moreover, we see that the d-wave phase is stabilized close 
to quarter filling, in contrast with the 2D case. 

In conclusion, we have studied the eflFect of e-ph cou- 
pling on the phase diagram of a two-leg ladder, using 
the Holstein-Hubbard model and applying the functional 
renormalization group method. This method deals with 
the retardation effects in a systematic way. In partic- 
ular, when the e-ph interaction strength is comparable 
to the e-e interactions and A ^ 1, there are qualitative 
differences in the phase diagram from that predicted in 
previous studies. In particular, in addition to the SDW 
and d-wave SC phases found previously, we establish that 
for strong e-ph coupling, both CDW and s-wave SC can 
arise. 

A. H.C.N, was supported through NSF grant DMR- 
0343790. 
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